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We study theoretically two-dimensional single-crystalline sheets of semiconductors forming a hon¬ 
eycomb lattice with a period below 10 nm. These systems could combine the usual semiconductor 
properties with Dirac bands. Using atomistic tight-binding calculations, we show that both the 
atomic lattice and the overall geometry influence the band structure, revealing materials with un¬ 
usual electronic properties. In rock-salt Pb-chalcogenides, the expected Dirac-type features are 
clouded by a complex band structure. However, in the case of zinc-blende Cd-chalcogenide semicon¬ 
ductors, the honeycomb nano-geometry leads to rich band structures including, in the conduction 
band, Dirac cones at two distinct energies and non-trivial flat bands, and, in the valence band, 
topological edge states. These edge states are present in several electronic gaps opened in the va¬ 
lence band by the spin-orbit coupling and the quantum confinement in the honeycomb geometry. 

The lowest Dirac conduction band has ^-orbital character and is equivalent to the 7r — 7r* band of 
graphene but with renormalized couplings. The conduction bands higher in energy have no counter¬ 
part in graphene, they combine a Dirac cone and flat bands because of their P-orbital character. We 
show that the width of the Dirac bands varies between tens and hundreds of meV. These systems 
emerge as remarkable platforms for studying complex electronic phases starting from conventional 
semiconductors. Recent advancements in colloidal chemistry indicate that these materials can be 
synthesized from semiconductor nanocrystals. 

I. INTRODUCTION 

The interest in two-dimensional (2D) systems with 
a honeycomb lattice and related Dirac-type electronic 
bands has exceeded the prototype graphene [1]. Cur¬ 
rently, 2D atomic @-[6| and nanoscale [7Hlljl systems are 
extensively investigated in the search for materials with 
novel electronic properties that can be tailored by geome¬ 
try. For example, a confining potential energy array with 
honeycomb geometry was created on a Cu(lll) surface 
and it was demonstrated that the electrons of the Cu 
surface state have properties similar to those of graphene 
[6]. From the same perspective, it was proposed that a 
honeycomb pattern with a 50-100 nm periodicity could 
be imposed on a 2D electron gas at the surface of a con¬ 
ventional semiconductor by using lithography or arrays 
of metallic gates Within the effective-mass ap¬ 

proach, linear E( k) relationships were predicted close to 
the Dirac points in the Brillouin zone, in analogy with 
graphene Q . The group velocity of the carriers was found 
to be inversely proportional to the honeycomb period and 
to the effective carrier mass. In order to obtain a system 
with sufficiently broad Dirac bands, it is thus of major 
importance to reduce the period of the honeycomb lat¬ 
tices far below 50 nm and to use semiconductors with low 
effective mass. 
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Our aim in the present work is to explore theoretically 
the physics of 2D semiconductors with honeycomb geom¬ 
etry and period below 10 nm. Electronic structure calcu¬ 
lations using the atomistic tight-binding method Bill 
attest that both the atomic lattice and the overall ge¬ 
ometry influence the band structure. We show that the 
honeycomb nano-geometry not only enables the realiza¬ 
tion of artificial graphene with tunable properties but 
also reveals systems with non-trivial electronic structure 
which has no counterpart in real graphene. 

We consider atomically-coherent honeycomb superlat¬ 
tices of rock-salt (PbSe) and zinc-blende (CdSe) semi¬ 
conductors. These artificial systems combine Dirac-type 
electronic bands with the beneficial tunability of semi¬ 
conductors under strong quantum confinement. In the 
case of a zinc-blende atomic lattice, separated conduc¬ 
tion IS and IP Dirac cones of considerable bandwidth 
(10’s to 100’s of meV) are found, as well as dispersionless 
IP bands. Here, IS and IP refer to the symmetry of 
the wave-functions on each node of the honeycomb. The 
chirality of the wave-functions with respect to apseudo- 
spin is also demonstrated for both Dirac cones [1]. This 
rich electronic structure is attributed to the absence of 
hybridisation between IS and IP bands. We show that 
the physics for fermions in honeycomb optical lattices of 
cold atoms with p orbitals [lj, H5| could be studied in 
nanostructured 2D semiconductors. We point out subtle 
differences between the electronic structure of graphene- 
and silicene-type honeycomb structures. In the latter 
case, gaps at the Dirac points can be controllably opened 
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FIG. 1. (a,b) Block models for the self-assembled honeycomb 
lattices based on nanocrystals that have a truncated cubic 
shape, (a) Honeycomb lattice formed by atomic attachment 
of 3 of the 12 {110} facets (light green/orange) corresponding 
to a graphene-type honeycomb structure (both sub-lattices in 
one plane). (b) Honeycomb lattice formed by attachment of 3 
of the 8 {111} facets (dark green/red), leading to the silicene- 
type configuration (each sub-lattice in a different plane). In 
both models, a (111) direction is perpendicular to the hon¬ 
eycomb plane as is experimentally observed. The arrow indi¬ 
cates the superlattice parameter d. (c) Top view of the unit 
cell of a graphene-type honeycomb lattice of PbSe nanocrys¬ 
tals (Pb atoms: grey; Se atoms: yellow). 


and closed by an electric field applied perpendicularly to 
the 2D structure. In the valence band of CdSe sheets, we 
demonstrate the existence of topological edge states in 
the electronic gaps opened by the spin-orbit coupling and 
the quantum confinement, supporting the recent work of 
Sushkov et al. using envelope-function theory fl6j |. Our 
atomistic calculations even predict multiple gaps with 
edge states. 


II. GEOMETRY OF THE 2D LATTICES 

The 2D crystals that we consider in the present work 
are inspired by recent experiments [13 discussed in 
Sec.Kl These experiments consist in the synthesis in a 
2D reactor plane of honeycomb sheets of PbSe by self- 
assembly and atomic attachment of (nearly) monodis- 
perse PbSe colloidal nanocrystals with a truncated cu¬ 
bic shape [i3 El. Due to facet-specific atomic bonding, 
atomically-coherent lattices with a honeycomb geometry 
and long-range periodicity are formed. Rock-salt PbSe 
lattices are transformed into zinc-blende CdSe lattices 
by cation exchange chemistry E3E1. In both cases, the 
(111) axis of the atomic lattice is perpendicular to the 
plane of the honeycomb sheet. 

The single-crystalline sheets that we consider are made 
of nano crystals arranged in honeycomb structure, as 
shown in Fig. [lj They can be seen as triangular lat¬ 
tices with a basis of two nanocrystals per unit cell form¬ 
ing by periodicity sub-lattices A and B. There are only 
two ways to assemble truncated nanocubes with a (111) 
body diagonal upright into a honeycomb lattice. The 
first one is to use three {110} facets perpendicular to 
the superlattice plane with angles of 120 degrees, result¬ 
ing in the honeycomb lattices presented in Fig. [l}i. In 
this structure all nanocrystal units are organized in one 
plane, i.e. equivalent to a graphene-type honeycomb lat¬ 
tice. In the second configuration, three truncated {111} 
planes per nanocrystal are used for atomic contact. This 
results in atomically crystalline structures in which the 
nanocrystals of the A and B sub-lattices are centred on 
different parallel planes, in analogy with atomic silicene 
(Fig. [TJd). The separation between the two planes is 
d/( 2\/6) ~ 0.2 d where d is the superlattice parameter 
defined in Fig. [TJ). The unit cell of a typical graphene- 
type honeycomb lattice of PbSe nanocrystals is displayed 
in Fig. |Tt, in the case where polar {111} facets are termi¬ 
nated by Pb atoms. We have also considered nanocrys¬ 
tals with Se-terminated {111} facets. Further details on 
the lattice geometry and nanocrystal shape related to the 
degree of truncation are given in the Appendix [A] 

The band structures presented below have been calcu¬ 
lated for these quite specific geometries inspired by ex¬ 
periments. However, very similar results can be obtained 
for other honeycomb nano-geometries, as shown in Ap¬ 
pendix o for honeycomb lattices composed of spherical 
nanocrystals connected by cylindrical bridges. 

III. TIGHT-BINDING METHODOLOGY 

Based on the effective-mass approach, Dirac-type 
bands of considerable width can be expected in a 
nanoscale honeycomb lattice jsj. However, a detailed un¬ 
derstanding of the electronic structure in relation to the 
atomic structure and nanoscale geometry of these sys¬ 
tems requires more advanced calculation methods. We 
have, therefore, calculated the energy bands of honey- 
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comb lattices using an atomistic tight-binding method. 
Each atom in the lattice is described by a double set 
of sp 3 d 5 s* atomic orbitals including the spin degree of 
freedom. We include spin-orbit interaction and we use 
tight-binding parameterizations (Appendix [Bj) that give 
very accurate band structures for bulk PbSe and CdSe. 
To avoid surface states, CdSe structures are saturated by 
pseudo-hydrogen atoms. This is not necessary in rock- 
salt PbSe structures as discussed in Ref. [20]. Due to 
the large size of the systems that we have studied (up 
to 6 x 10 4 atoms and 1.2 x 10 6 atomic orbitals per unit 
cell), only near-gap eigenstates are calculated using the 
numerical methods described by Niquet et al 0. 


IV. CONDUCTION BAND STRUCTURE OF 
GRAPHENE-TYPE LATTICES OF CdSe 

The typical dispersion [P(k)] of the highest occupied 
bands and the lowest unoccupied bands of a graphene- 
type superlattice of CdSe is shown in Fig. [2] The elec¬ 
tronic structure is composed of a succession of bands and 
gaps due to the nanoscale periodicity. The honeycomb 
geometry induces periodic scattering of the electronic 
waves, opening gaps in particular at the center and at 
the edges of the superlattice Brillouin zone. 

We first discuss the simpler but extremely rich physics 
of the lowest conduction bands (Fig.[2ji), which consist of 
two well-separated manifolds of two and six bands (four 
and twelve bands including spin). Strikingly, the two low¬ 
est bands have the same type of dispersion as the tt and 
7r* bands in real graphene; these bands are connected 
just at the K (k = K) and K’ points of the Brillouin 
zone, where the dispersion is linear (Dirac points). More¬ 
over, in the second manifold higher in energy, four bands 
have a small dispersion and two others form very disper¬ 
sive Dirac bands. The presence of two separated Dirac 
bands with Dirac points which can be experimentally ac¬ 
cessed by the Fermi level is remarkable and has not been 
found in any other solid-state system. It can be under¬ 
stood from the electronic structure of individual CdSe 
nanocrystals characterized by a spin-degenerate electron 
state with a IS envelope wave-function and by three spin- 
degenerate IP excited states higher in energy. The two 
manifolds of bands arise from the inter-nanocrystal cou¬ 
pling between these IS and IP states, respectively. In¬ 
terestingly, with this nanoscale geometry, the coupling 
between nanocrystal wave-functions is strong enough to 
form dispersive bands with high velocity at the Dirac 
points, but small enough to avoid mixing (hybridization) 
between IS and IP states. 

The systematic presence of nearly flat IP bands is an¬ 
other remarkable consequence of the absence of S — P 
hybridization. The existence of dispersionless bands has 
been predicted in honeycomb optical lattices of cold 
atoms with p orbitals [lj, [l5j. I n our case, two IP bands 
are built from the 1 P z states perpendicular to the lattice, 
they are not very dispersive simply because 1P Z -1P Z (i r) 



FIG. 2. Lowest conduction bands (a) and highest valence 
bands (b) of a graphene-type honeycomb lattice of truncated 
nanocubes of CdSe (body diagonal of 4.30 nm). The zero of 
energy corresponds to the top of the valence band of bulk 
CdSe. (c) Energy splittings (meV) between the two highest 
valence bands (VB), between the two IS Dirac conduction 
bands (IS), and between the two IP Dirac conduction bands 
(IP) calculated at K + q and plotted using polar coordinates, 
i.e., energy versus the angle between K and q (squares, red 
curve: |q| — 0.05|K| ; circles, blue curve: |q| as 0.15|K| ; 
triangles, black curve: |q| = 0.25|K|). 

interactions are weak. Two other IP bands (1P Xj2/ ), re¬ 
spectively above and below the IP Dirac band, are flat 
due to destructive interferences of electron hopping in¬ 
duced by the honeycomb geometry mint (non-trivial 
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FIG. 3. Evolution of the bandwidth of the IS Dirac band as 
a function of the number of atoms that form the nanocrys¬ 
tal/nanocrystal contact (./Vat ) • The size of the nanocrystals, 
which also determines the period of the honeycomb lattice, 
is indicated by the body diagonal of the truncated nanocube. 
Results are shown for the graphene geometry (bonding via 
the {110} facets). 


flat bands). 

Close to the Dirac points, for k = K + q where 
|q| < 0.1 |K|, the dispersion of the IS and IP Dirac bands 
is remarkably isotropic, i.e., it does not depend on the an¬ 
gle between q and K (Fig. [2k). For larger values of |q|, 
the Dirac cones exhibit trigonal warping due to the effect 
of the super lattice potential on the electrons. Trigonal 
deformation of the energy bands has a profound impact 
on the (quantum) Hall effect, interference patterns and 
weak localization in graphite [22] and bilayer graphene 
[23j . Unusual phenomena such as enhanced interference 
around defects and magnetically ordered exotic surfaces 
are also predicted at the surface of 3D topological insula¬ 
tors due to the hexagonal warping of the bands mm. 
By analogy, similar effects could arise in the honeycomb 
super lattices of nanocrystals. 

At T (k = 0), the width of the IS band increases 
not only with decreasing nanocrystal size but also with 
increasing number of contact atoms (Fig. [3]). This can be 
understood by the fact that the contact area determines 
the electronic coupling between the nanocrystal wave- 
functions of adjacent sites, and acts in a similar way as 
the hopping parameter in atomic honeycomb lattices (see 
Sec. US]). We predict a bandwidth above 100 meV for 
realistic configurations. Note that the width of the Dirac 
IP bands is even considerably larger than that of the 
IS bands. In the case of the graphene geometry, we have 
found that the IS and IP bands are always characterized 
by well-defined Dirac points. 

Conduction bands at higher energy (> 2.6 eV) in 
Fig. [2k are derived from the ID wave-functions of the 
nanocrystals. Interestingly, they also present flat bands 
induced by the honeycomb geometry. 


V. NON-TRIVIAL GAPS IN THE VALENCE 

BAND OF GRAPHENE-TYPE LATTICES OF 
CdSe 

The twofold-degenerate conduction band of bulk CdSe 
near the T point is mainly derived from s atomic orbitals 
and is therefore characterized by a very weak spin-orbit 
coupling. It is the reason why IP and IP bands of CdSe 
superlattices are spin-degenerate and exhibit well-defined 
Dirac points. The situation is totally different in the va¬ 
lence bands of CdSe which are built from p 3 / 2 — P 1/2 
atomic orbitals characterized by a strong spin-orbit cou¬ 
pling, leading in the bulk to a splitting of 0.39 eV between 
heavy-hole and split-off bands at T. 

Therefore graphene-type lattices of CdSe nanocrystals 
present two extremely interesting features when they are 
combined, 1) a honeycomb geometry, 2) a strong spin- 
orbit coupling, here in the valence band. In this kind of 
systems, the spin-orbit coupling may open a non-trivial 
gap and give rise to topological insulators. These sys¬ 
tems exhibit remarkable properties at their boundaries 
characterized by helical edge states in the gap induced 
by the spin-orbit coupling [26l428j| . The edge states carry 
dissipation-less currents, leading in 2D systems to quan¬ 
tum spin Hall effect, as initially predicted for graphene 
by Kane and Mele [29|. Whereas in graphene the spin- 
orbit coupling is too small to give measurable effects [1], 
it was predicted that topological insulators with much 
larger gaps could be made from ordinary semiconductors 
on which ^potential with hexagonal symmetry is super¬ 
imposed [16|. We show below that honeycomb lattices 
of CdSe nanocrystals actually present several non-trivial 
gaps in their valence band. 

The dispersion of the valence bands in CdSe super¬ 
lattices (Fig. [ 2 ) 3 ) is much more complex than the disper¬ 
sion of conduction bands because anisotropic heavy-hole, 
light-hole and split-off bands are coupled by the confine¬ 
ment. Dirac points are not visible in spite of the honey¬ 
comb geometry. However, the highest valence bands in 
Fig. [ 2)3 roughly behave like the tt — tt* bands in graphene, 
their energy-momentum relationship is isotropic close to 
the K point (Fig. [2k), but there is a large gap between 
the two bands. Other gaps are present at lower energy. 

In the following, we investigate the topological proper¬ 
ties of these gaps. First, we calculate the Chern numbers 
for the two highest valence bands of the 2D sheet (30j . 
The methodology and the results are described in Ap¬ 
pendix [Dj The analysis of these Chern numbers demon¬ 
strates the non-trivial character of the two gaps between 
the highest valence bands and therefore the quantum spin 
Hall effect is predicted in these gaps. Second, we calcu¬ 
late the valence band structure of ribbons with zigzag 
and armchair geometry built from CdSe nanocrystals. 
We have considered the same honeycomb geometry as in 
Fig. [ 2)3 but for a ID ribbon instead of a 2D sheet. The 
unit cell forming the ribbon by periodicity is composed 
of sixteen nanocrystals. Figure [4k presents the ID band 
structure for a ribbon with zigzag edges. We have found 
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FIG. 4. Highest valence bands in ribbons made from a 
graphene-type honeycomb lattice of truncated nanocubes of 
CdSe (body diagonal of 4.30 nm). (a) Ribbon with zigzag 
edges (ribbon width = 57 nm, periodic cell length l — 8.2 
nm). (b) Ribbon with armchair edges (ribbon width = 33 
nm, periodic cell length l — 14.2 nm). In each case, the unit 
cell is composed of sixteen nanocrystals (34768 atoms per unit 
cell). The coloured regions indicate the bands of the corre¬ 
sponding 2D semiconductor, i.e., those shown in Fig. [2)3. 
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FIG. 5. 2D plots of the wave-functions of the edge states 
calculated at k = 0.3x27r/Z for the ribbon considered in FigHJi 
(energy of the states ~ —0.123 eV). The plots are restricted 
to a single unit cell of the ribbon (a: higher energy state; 
b: lower energy state). The vertical axis corresponds to the 
direction of the ribbon. More than 90% of each wave-function 
is localized on the nanocrystals at the edges. The white dots 
indicate the atoms. 



FIG. 6. Evolution of the energy gap between the two highest 
valence bands of CdSe sheets with graphene-like honeycomb 
geometry as a function of the number of atoms that form 
the nanocrystal/nanocrystal contact (JV a t). The size of the 
nanocrystals is indicated by the body diagonal of the trun¬ 
cated nanocube. 


edge states crossing the gap between the two highest va¬ 
lence bands of the 2D sheet. The 2D plot of their wave- 
functions is shown in Fig. [5] Interestingly, edge states 
are also present in the second gap below in energy, and 
others are even visible in the smaller gaps below -0.15 eV. 
The analysis of their wave-functions on one side of the 
ribbon shows that the spin is mainly oriented perpendic¬ 
ular to the lattice (> 98%) [31], and that the direction 
of the spin is reversed for motion in the opposite direc¬ 
tion (fc —)► —k). The situation is inverted at the opposite 
edge of the ribbon. Very similar results are obtained for 
ribbons with armchair edges (Fig. [4 )d) and for ribbons in 
which we modify the size of the nanocrystals at the edges 
(not shown). The presence of helical edge states in the 
gaps of the 2D sheet and their robustness with respect 
to the edge geometry are signatures of their non-trivial 
topology. The multiplicity of non-trivial gaps demon¬ 
strates the variety of effects induced by the honeycomb 
nano-geometry on the band structure of the 2D semicon¬ 
ductor, contrarily to the case of a square nano-geometry 


The gaps between the valence bands are tunable 
thanks to the quantum confinement. Figure [6] shows that 
the gap between the two highest valence bands strongly 
depends on the size of the nanocrystals and increases 
with the number of atoms at the contact plane between 
neighbor nanocrystals. A gap above 10 meV is possible 
with nanocrystal size below 4 nm. 


VI. PSEUDO-SPIN 

A fingerprint of the electronic states at the Dirac cones 
in graphene is their chirality with respect to a pseudo¬ 
spin associated with the two components of the wave- 
function on the two atoms of the unit cell [il]. Figure 0 
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FIG. 7. Chirality of the IP wave-functions of the CdSe honey¬ 
comb lattice compared with atomic chirality in graphene, (a) 
Phase shift (radians) on each atomic orbital between the elec¬ 
tronic states of the lowest conduction band of the graphene- 
type honeycomb superlattice of the compound CdSe (Fig.[2jt). 
The phase shift is calculated at K + qi and K + q 2 with 
|qi| = |q 2 1 as 0.05|K| and an angle between qi and q 2 of 3.8 
radians, for each lateral atomic position in the lattice plane, 
(b) Difference between the phase of the wave-function at the 
center of nanocrystals A and B versus the angle between qi 
and q 2 , at K (blue, square symbols) and K’ (red, circular 
symbols). Similar results are obtained for the IP Dirac cone. 


shows that the pseudo-spin is also well defined near the 
Dirac points in the conduction band of CdSe honeycomb 
lattices, in spite of the fact that each unit cell of the 
structure contains thousands of atoms. When we rotate 
the k vector around the K point, the phase shift of the 
wave-function is almost constant across each nanocrystal, 
in other words, it does not depend on the atomic orbital 
and its position. Note that the wave-function phase shift 
changes quite abruptly at the contact plane between two 
nanocrystals (one of sub-lattice A and one of sub-lattice 
B). In the lower IS band, the phase difference between 
nanocrystals A and B is equal to the angle of rotation, 
and the variation is opposite at the K’ point (Fig. [7t>). 
The sign is also inverted in the upper IS band. This chi¬ 
rality of the wave-function should reduce backscattering 
of the Dirac electrons, for the same reasons as in graphene 

Q. 



Electric field (meV/nm) 

FIG. 8. (a) Lowest conduction bands of a silicene-type honey¬ 
comb lattice of truncated nanocubes of CdSe (body diagonal 
of 5.27 nm). (b) Evolution of the gap at K between the IS 
(blue squares) and between the IP (red circles) Dirac bands 
versus the electric field strength applied perpendicular to the 
honeycomb plane. 


VII. BAND STRUCTURE OF SILICENE-TYPE 
LATTICES OF CdSe 


The band structures for honeycomb lattices of CdSe 
with silicene-type geometry are almost the same as for 
the graphene-type geometry (see Fig. [8fct for conduction 
bands), but there are gaps at the Dirac points due to the 
absence of mirror symmetry with respect to the {111} 
contact plane between neighbour nanocrystals. In other 
words, sub-lattice A is not equivalent to sub-lattice B. 
The gap at K strongly depends on the size and the shape 
of the nanocrystals. Since in that case the nanocrys¬ 
tals of the A and B sub-lattices are positioned at differ¬ 
ent heights, the electronic structure is very sensitive to 
an electric field applied along (111) as shown in Fig [HJ d. 
The gap at K between the IP bands varies linearly with 
the field and vanishes when the potential drop between 
nanocrystals A and B compensates the effect of the ge¬ 
ometrical asymmetry on the IP wave-functions. The 
variation of the gap at K between IP Dirac bands is 
more complex due to intermixing between IP nanocrys¬ 
tal wave-functions, and increasing IP — IP and IP — ID 
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FIG. 9. Lowest conduction bands (a) and highest valence 
bands (b) of a graphene-type honeycomb lattice of truncated 
nanocubes of PbSe (body diagonal of 5.00 nm, q = 0.5). 

hybridization. However, this gap tends to zero at increas¬ 
ing electric field. 

VIII. BAND STRUCTURE OF 
GRAPHENE-TYPE LATTICES OF PbSe 

Figure [9] shows a typical band structure of a honey¬ 
comb lattice of PbSe with graphene-type geometry. The 
lowest conduction bands and highest valence bands are 
characterized by a manifold of eight bands which are 
formed by the 2 x 4 IS conduction and valence wave- 
functions of the two PbSe nanocrystals of the unit cell. 
These IS nanocrystal wave-functions are derived from 
the fourfold degenerate conduction and valence bands of 
bulk PbSe at the L point of the Brillouin zone (eightfold 
degeneracy if we include the spin). The conduction bands 
at higher energy and the valence bands at lower energy 
are derived from the 24 IP nanocrystal wave-functions 
(only seven bands are shown for clarity). The complex 
dispersion of all these bands, even the IS ones, shows 
that the coupling between the wave-functions of neigh¬ 
bor nanocrystals is not only governed by the symmetry 
of the envelope function (1 S', IP) defined by the nano¬ 
geometry but also, in a subtle way, by the underlying 
Bloch function which depends on the originating valley. 


These results illustrate the fact that both the atomic lat¬ 
tice of the parent semiconductor and the overall nano¬ 
geometry influence the band structure. 

IX. EFFECTIVE TIGHT-BINDING MODEL 

The overall behavior of the IP and IP conduction 
bands of honeycomb lattices of CdSe nanocrystals can 
be interpreted using a simple (effective) tight-binding 
Hamiltonian of graphene or silicene in which the two 
” atoms” of the unit cell are described by one 5 and three 
p orbitals (doubled when the spin is considered). The on¬ 
site energies are E S (A), E Px (A ) = E Py (A ), E Pz (A ), and 
E S (B ), E Px (B ) = E Py (B ), E Pz (B) for the A and B sub¬ 
lattices, respectively. The 2 axis is taken perpendicular 
to the lattice. For silicene, sub-lattice B is not in the same 
plane as A, the angle between the AB bond and the z axis 
is taken to be the angle between the [111] and [11 — 1] crys¬ 
tallographic axes. The energy of the p z orbital is allowed 
to be different from the p x and p y orbitals. Following 
Slater and Koster [33j, all nearest-neighbor interactions 
(hopping terms) can be written in the two-center approx¬ 
imation as functions of four parameters (V SS(T , V spa , V ppa , 
V pp7r ) plus geometrical factors. The problems for the s 
and p bands are separable when V spa = 0. We have found 
that the results of full tight-binding calculations for all 
2D crystals considered are only compatible with small 
values of V spcri i.e., by a small hybridization between s 
and p orbitals. 

Typical band structures for the silicene geometry are 
shown in Fig [lOl Figure [TOb corresponds to a case where 
the on-site energies are the same on the sub-lattices A 
and B, and all p orbitals have the same energy. Dirac 
points are obtained at three energies, one in the 8 bands, 
two in the p, and no gap formed at K. This situation has 
not been found for more realistic band structure calcula¬ 
tions. One Dirac p band is suppressed when we slightly 
shift the p z energy with respect to the p x , p y energies 
(Fig flQb). Finally, a small gap is opened at the two 
Dirac points when the orbital energies in sub-lattice A 
are not the same as in B (FigfTOb). Figure [TOb is in close 
agreement with Fig. [8k highlighting the suitability of the 
present effective Hamiltonian to describe the obtained 
results. 


X. EXPERIMENTAL PERSPECTIVE 

Several methods can be proposed to confine semicon¬ 
ductor carriers in a honeycomb geometry. The simplest 
method would be to use a hexagonal or honeycomb array 
of a metal acting as a geometric gate that would electro¬ 
statically force the carriers in a honeycomb lattice. The¬ 
oretical and experimental work with GaAs/GaAlAs has 
been recently reported [9, [I(J [n| . The advantage of this 
approach is the substantial theoretical and experimental 
knowledge already attained in the field. However, due to 
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FIG. 10. Band structure calculated for a silicene-type lattice 
using the effective tight-binding Hamiltonian for three sets of 
parameters, (a) E S (A : B) = 2.18 eV, all the E p = 2.44 eV, 
Vsscr — 8 meV, Vspa — 0 meV, V PP a — 50 meV, V PP n — 2 

meV. (b) Same as (a) but with E Pz (A,B ) = 2.51 eV. (c) 
E S (A) = 2.18 eV, E S (B) = 2.19 eV, E Px (A) = E Py (A) = 
2.44 eV, E Px (B) = E Py (B ) = 2.46 eV, E Pz (A) = 2.51 eV, 
E Pz (B) — 2.52 eV, V SS a = -8 meV, V spa = 0 meV, V ppa — 50 
meV, V pp7T = — 2 meV. 


the relatively large period of the honeycomb lattice, the 
width of the bands is limited to a few meV [H. A second 
method would be to prepare suitable templates by lithog¬ 
raphy, and use these to grow semiconductor lattices with 
wet-chemical or gas-phase methods (chemical vapor de¬ 
position and molecular beam epitaxy). General concerns 
relate to the crystallinity of the 2D system prepared with 
a honeycomb geometry. Our tight-binding calculations 
unambiguously show that the electronic structure is de¬ 
termined both by the periodicity of the 2D semiconduc¬ 
tor in the nanometer range and the atomic lattice. Re¬ 
cently, a third approach to prepare semiconductors with 


a honeycomb geometry came from a rather unexpected 
corner. Evers et al. showed that atomically-coherent 
sheets of a PbSe semiconductor can be prepared by ori¬ 
ented attachment of colloidal nanocubes [13]. Moreover, 
the as-prepared systems can be transformed into zinc- 
blende CdSe by Cd-for-Pb ion exchange, preserving the 
nano-geometry [34j] . The systems showed astounding ge¬ 
ometrical order and a well-defined atomic structure. A 
challenge with theses systems will be to incorporate them 
in a field-effect transistor, such that the transport prop¬ 
erties can be measured. As with the template-approach, 
the effect of defects in the atomic lattice and disorder in 
the superimposed honeycomb geometry, electronic dop¬ 
ing, and the surface termination of the 2D honeycomb 
semiconductors are issues that need to be further ad¬ 
dressed. 


XI. SUMMARY AND FUTURE DIRECTIONS 


In conclusion, we have shown that atomically-coherent 
honeycomb lattices of semiconductors with period below 
10 nm present very interesting band structures, which 
are defined both by the properties of the parent semi¬ 
conductor and the nano-geometry. We have calculated 
the electronic structure of these materials using an atom¬ 
istic tight-binding method, considering graphene- and 
silicene-type lattices obtained by the assembling of semi¬ 
conductor nanocrystals. In the case of honeycomb lat¬ 
tices of CdSe, we predict a rich conduction band struc¬ 
ture exhibiting non-trivial flat bands and Dirac cones 
at two experimentally reachable energies. These bands 
are derived from the coupling between nearest-neighbor 
nanocrystal wave-functions with IS and IP symmetry. 
The formation of distinct Dirac cones is possible be¬ 
cause of the weak hybridization between IS and IP wave- 
functions under the effect of the strong quantum confine¬ 
ment. We also predict the opening of non-trivial gaps in 
the valence band of CdSe sheets due to the effect of the 
nano-geometry and the spin-orbit coupling. Several topo¬ 
logical edge states are found in these gaps. The possibil¬ 
ity to have multiple Dirac cones, non-trivial flat bands 
and topological insulating gaps in the same system is re¬ 
markable. 

Recent experiments strongly suggest that the synthe¬ 
sis of such single-crystalline sheets of semiconductors is 
possible using facet-specific attachement of nanocrystals 
(17j |. Therefore, numerous directions are open for the 
theoretical and experimental investigation of these sys¬ 
tems. The electronic structure and carrier transport can 
be studied by local scanning tunnelling microscopy and 
spectroscopy, and in a field-effect transistor. With illu¬ 
mination and/or gating the conduction band can be con- 
trollably filled with electrons up to several electrons per 
nanocrystal [3514371] , allowing for the Fermi level to cross 
the Dirac points. Our work provides evidence for non¬ 
trivial flat IP bands, and the lowest one can be reached at 
a nanocrystal filling between 2 and 3 electrons. Electron- 
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electron interactions should then play a crucial role and 
may lead to Wigner crystallization [hj] . We remark that 
the physics of honeycomb lattices of p orbitals is largely 
unexplored. It will be particularly attractive to consider 
honeycomb lattices of semiconductors with even stronger 
spin-orbit coupling, in which we expect the emergence of 
new electronic (topological) phases. 

Finally, it is important to realize that all these in¬ 
teresting properties (Dirac cones, non-trivial flat bands 
and topological edge states) originate from the honey¬ 
comb geometry. For instance, the band structures of 
atomically-coherent square semiconductor superlattices 
that we discussed recently [32], although of interest in 
their own respect, do not contain any of the Dirac-based 
quantum electronic properties of similar semiconductors 
with a honeycomb nano-geometry, discussed here. 
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Appendix A: Details on the geometry of the 
honeycomb lattices 

1. Formation of the nanocrystals 

The nanocrystals are built from PbSe nanocubes on 
which six {100}, eight {111} and twelve {110} facets are 
created by truncation. The positions of the vertices of the 
nanocrystal shape are given by P[d=l, ±(1 — g), ±(1 — q)\ 
where [d=l, =bl, ±1] indicate the position of the six corners 
of the original nanocube, q is the truncation factor, and 
P represents all the possible permutations. We have con¬ 
sidered realistic shapes corresponding to q between 0.25 
and 0.5 E3- The truncated (lll)-oriented nanocubes 
are assembled into two types of honeycomb lattices with 
a graphene-like or silicene-like shape as described below. 


2. Geometry of graphene-type supercells 

Two nanocrystals (A and B) are attached along the 
(110) direction (perpendicular to the (111) axis) to define 
a unit cell, forming by periodicity an atomically-coherent 
honeycomb lattice. Each nanocrystal is defined by the 
same number of bi-planes of atoms in the (110) direction 
but there is an additional plane of atoms shared between 
neighboring nanocrystals in order to avoid the formation 
of wrong bonds. The length of the two vectors defining 
the superlattice is (2n + l)a>/(6)/2 where n is an integer 


TABLE I. Tight-binding parameters (notations of Slater and 
Foster Hi) for zinc-blende CdSe in an orthogonal sp 3 d 5 s* 
model. A is the spin-orbit coupling, (a) and (c) denote the 
anion (Se) and the cation (Cd), respectively. 


Parameters for CdSe (eV) 


E s (a) 

-8.065657 

Es( c) 

-1.857148 

Ep (a) 

4.870028 

E p (c) 

5.613460 

Ed xy (a) 

15.671502 

Ed xy ( c) 

16.715749 

E d x 2_ y 2 ( a ) 

15.232107 

E d x 2_ y 2 ( C ) 

20.151047 

E s * (a) 

15.636238 

E s *( c) 

20.004452 

A(a) 

0.140000 

A(c) 

0.150000 

Esso- (ac) 

-1.639722 

Us* s* cr (ac) 

-1.805116 

Uss* a (ac) 

1.317093 

Uss*CT (ca) 

0.039842 

Uspcr (ac) 

3.668731 

Uspcr (ca) 

1.885956 

Us*per (ac) 

0.978722 

Vs*per (ca) 

1.424094 

Usdcr (ac) 

-0.890315 

Usdcr (ca) 

-1.007270 

Us* da (ac) 

0.906630 

Us* dcr (ca) 

2.472941 

Upper (ac) 

4.430196 

Upp 7 r (ac) 

-0.798156 

Vpdcr (ac) 

-2.645560 

Vpdcr (ca) 

-1.296749 

Upd 7 r (ac) 

0.028089 

Upd 7 r (ca) 

2.295717 

Udder (ac) 

-2.480060 

Udd 7 r (ac) 

2.393224 

Vdd5 

-1.373199 




Parameters for Cd-H and Se-H (eV) 

E h 0.000000 

V 8S <t -35.69727 V sp(T 61.82948 


and a is the lattice parameter (0.612 nm for PbSe, 0.608 
nm for CdSe). Since the atomistic reconstruction at the 
contact plane between neighboring nanocrystals is not 
precisely known, we have also considered structures in 
which we have slightly enlarged each inter-nanocrystal 
junction by one line of atoms on each side of the {110} 
facets. 


3. Geometry of silicene-type supercells 

The two nano crystals A and B are attached along 
a {111} facet to define a unit cell and an atomically- 
coherent honeycomb lattice is formed by periodicity. 
Note however that the A and B nanocrystal sub-lattices 
are each located in a different plane. The length of the su¬ 
perlattice vectors (superlattice parameter) is d = mna/y/2 
where m is an integer. The separation between the two 
planes A and B is raa/(4\/3). Similar to the graphene- 
type lattices, we have also considered structures in which 
the {111} facets at the contact plane between neighbor¬ 
ing nanocrystals are enlarged. 
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FIG. 11. Lowest conduction bands of a graphene-type hon¬ 
eycomb lattice of spherical nanocrystals of CdSe (diameter 
= 4.7 nm). Nearest-neighbor nanocrystals are connected by 
cylindrical bridges of CdSe (diameter of the cylinders = 2.3 
nm). 


Appendix B: Tight-binding parameters 

We consider a double basis of sp 3 d 5 s* atomic orbitals 
for each Pb, Cd or Se atom, including the spin degree of 
freedom. For PbSe, we use the tight-binding parameters 
as given in Ref. [20|. Due to the lack of sp 3 d 5 s* tight- 
binding parameters for zinc-blende CdSe, corresponding 
data were derived and are presented in Table [H The use 
of a sp 3 d 5 s* basis allows us to get a very reliable band 
structure for bulk materials compared to ab initio cal¬ 
culations and available experimental data: see Ref. [20] 
for PbSe. Hence, these parameters can be safely trans¬ 
ferred to predict the electronic structure of semiconduc¬ 
tor nanostructures [13]. 

The effects of the electric field on the electronic struc¬ 
ture are calculated assuming that the field inside the 
superlattice is uniform. In this approximation, the cal¬ 
culation of the screening in such a complex system is 
avoided. The main conclusions of these calculations, i.e., 
the strong evolution of the bands with the applied field 
and the possibility to close or open the gap at the IS 
Dirac point in silicene-type lattices, do not depend on 
this approximation. 


Appendix C: Results for another type of honeycomb 
nano-geometry 

It is also possible to build a honeycomb lattice us¬ 
ing spherical nanocrystals. We set the distance between 
nearest-neighbor nanocrystals equal to the diameter, i.e., 
the spheres are tangential. Between each pair of neigh¬ 
bors, we add a cylinder of atoms which serves as a bond. 
We predict that the electronic bands for these 2D crys¬ 
talline sheets are very close to those obtained with trun¬ 
cated nanocubes. A typical result for the conduction 
band is shown in Fig. Qj] where the two Dirac cones and 
the flat bands are present. We have found similar band 


structures for diameters of the cylinders up to 80% of the 
diameter of the spheres, demonstrating the robustness of 
the results. The width of the Dirac bands is related to 
the strength of inter-nanocrystal bonds, i.e., to the di¬ 
ameter of the cylinders. The non-trivial flat bands are 
easily identified in Fig. |TT] since they are connected to 
the IP Dirac band. Similar results are also obtained for 
silicene-type geometry (not shown). 


Appendix D: Calculation of the Chern numbers 

We have calculated the Chern numbers of the high¬ 
est valence bands following the methodology proposed 
by Fukui et al. [38|. For each band n, we calculate the 
wave-functions |n, k) on a N x N grid within the Brillouin 
zone, i.e., k = raia*/7V + where a * and are 

the reciprocal lattice vectors (mi, m 2 = 0, ...,7V — 1). In 
order to specify the gauge, we make the transformation 
|n, k> |(n, k|0)| 1 (n, k|0)|n, k) where |0) is an arbi¬ 
trary state which was chosen as constant over the unit 
cell. The lattice Chern number associated with the nth 
band is calculated as 


< D1 > 

k 


where F( k) is defined as 


F(k) = In 


U!(k)U 2 




^2 _1 (k) 


(D2) 

with Vi (k) = (n, k|n, k + a*/7V)/|(n, k|n, k + a*/7V)|. 
F(k) is defined within the principal branch of the loga¬ 
rithm (— 7 r < F(k)/i < 7 r). It was shown that the lattice 
Chern number c n tends toward the usual Chern number 
c n in the limit N 00 [38j . 

In the calculation of the Chern numbers for the super¬ 
lattices, numerical and fundamental problems may arise. 
Numerical problems come from the size of the systems 
containing thousands of atoms per cell. However, the 
calculation of the quantities P(k) can be done separately 
for each value of k on the grid, enabling powerful paral¬ 
lel treatment on multi-processor computers. As a conse¬ 
quence, the methodology proposed by Fukui et al. [38] 
was found to be extremely efficient. 

Fundamental issues are more problematic. Each band 
among the highest valence bands shown in Fig. [2]3 is ac¬ 
tually composed of two almost-degenerate bands which 
intersect at several points of the Brillouin zone (T, K, M 
in particular). Elsewhere in the Brillouin zone, the split¬ 
ting between the almost-degenerate bands is smaller than 
0.2 meV in average (not visible in Fig. [ 2 ) 3 ). Therefore the 
methodology to calculate c n cannot be applied to these 
situations as it requires non-degenerate bands over the 
full Brillouin zone j38[. In order to lift the degeneracies 
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FIG. 12. Highest valence bands of a graphene-type honey¬ 
comb lattice of truncated nanocubes of CdSe (body diagonal 
of 4.30 nm) under a magnetic field B which induces Zeeman 
splitting of the bands (/isgsB — 1 meV). The corresponding 
bands at zero magnetic field are those shown in Fig. [2 Jd. The 
Chern numbers of the bands are indicated on the figure. 


between the bands, we have studied the effect of a mag¬ 
netic field applied perpendicularly to the superlattices. 
Since our objective is to explore the topological proper¬ 
ties of the bands, we have only considered the Zeeman 
part of the coupling Hamiltonian. Figure [12] shows that 
the application of a magnetic field totally splits the two 
highest valence bands, pushing downward (upward) the 
states with a majority of spin up (down) component. 

The calculated Chern numbers for the two highest 
manifolds of valence bands are indicated in Fig. [12] For 
each band, the sum in Eq. m typically converges to 
its final value for N ~ 30. We have checked that c n re¬ 
mains constant for larger N from « 30 to 100. It also 
remains constant when we vary the magnetic field, from 


HBgsB ~ 0.2 meV (minimum value to remove the dege- 
naracies) to 5 meV. It was not possible to calculate the 
Chern numbers for the valence bands lower in energy in 
this manner, because the degeneracy points remain even 
under an applied magnetic field. An equivalent approach 
to characterize the Z 2 topological invariant is given in 
Refs. dHH. 

For each valence band shown in Fig.[12j the sum of the 
Chern numbers c^ n + c^ n for the two spin components 
vanishes meaning that their contribution to the Hall con¬ 
ductivity is zero. However, the difference c^ n — c^ ?n (spin 
Chern number) which is linked to the spin Hall conduc¬ 
tivity m does n °t vanish. Since the sum of the Chern 
numbers over all the occupied valence bands must be 
zero for each spin component (there is no edge state in 
the gap between valence and conduction bands), we de¬ 
duce by subtraction that the quantum spin Hall effect 
must be present in the two highest gaps of the valence 
band, in total agreement with our previous conclusions 
of Sec. [Vl 

A small remark is in place here. For the reason that 
the spins are not pointing exactly up or down, the dif¬ 
ference c^ n — c^ n is strictly speaking not well defined. 
However, the spins are tilted only slightly away from the 
vertical (deviation < 2%), so that they can be treated 
as approximately up or down. The Chern numbers indi¬ 
cate that the corresponding spin “up” and spin “down” 
edge states counter propagate, i.e., that the bulk gap ex¬ 
hibits a nonzero spin Hall conductivity while the charge 
Hall conductivity vanishes. The small tilt of the spins 
perturbs the value of the spin Hall conductivity slightly 
away from the quantized value, but does not destroy the 
topological character of the state. Thus, this state is a 
quantum spin Hall state in good approximation. 
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